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Abstract: In spite of the advantages associated with the molecular 
specificity of fluorescence imaging, there is still a significant need to 
augment these approaches with label-free imaging. Therefore, we have 
implemented a form of interference microscopy based upon phase-shifted, 
laser- feedback interferometry and developed an algorithm that can be used 
to separate the contribution of the elastically scattered light by sub-cellular 
structures from the reflection at the coverslip-buffer interface. The method 
offers an opportunity to probe protein aggregation, index of refraction 
variations and structure. We measure the topography and reflection from 
calibration spheres and from stress fibers and adhesions in both fixed and 
motile cells. Unlike the data acquired with reflection interference contrast 
microscopy, where the reflection from adhesions can appear dark, our 
approach demonstrates that these regions have high reflectivity. The data 
acquired from fixed and live cells show the presence of a dense actin layer 
located 100 nm above the coverslip interface. Finally, the measured 
dynamics of filopodia and the lamella in a live cell supports retrograde flow 
as the dominate mechanism responsible for filopodia retraction. 
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1. Introduction 

Since its inception to measure the adhesion of cells to a transparent substratum, [1-3], inter- 
ference reflection microscopy (IRM) and reflection interference contrast microscopy (RICM) 
have yielded significant insights about the location and size of adhesion complexes and "focal 
contacts" (a name coined from an early application [4]). Although RICM can yield quantitative 
results when applied to imaging spheres [5] and vesicle adhesion [6, 7], the method, which is 
applicable to live cell imaging [8], has unfortunately not realized the potential for high precision 
measurements when applied to cells. One significant impediment is that analysis of the inten- 
sity variations in a set of interference fringes can lead to ambiguities because it is difficult to 
separate changes in topography from local variations in protein concentration near the ventral 
plasma membrane [9, 10]. Alternatively, the application of fluorescence interference contrast 
microscopy [9, 11, 12] has demonstrated that it is possible to use specific fluorescent label- 
ing to quantitatively probe the dynamics of membranes and the position of the ventral plasma 
membrane with nanometer precision. Indeed the molecular specificity of fluorescence based 
methods have enabled the imaging of the morphology of adhesions with resolution below the 
diffraction limit [13-15]. In spite of the advantages associated with fluorescent imaging, there 
is still a need to augment these approaches with label-free, interference imaging that can probe 
protein aggregation, index of refraction variations and structure [16, 17]. 

An approach based upon phase-shifting interferometry would mitigate these inherent limi- 
tations, but the potential impediments have appeared daunting because of the low reflectivity 
and complexity associated with elastic light scattering from sub-cellular features. Recent work, 
however, has demonstrated the power of using elastic light scattering methods [18] to probe 
actin dynamics [19,20], cell structure [21,22] and surface membrane waves [23,24] in live 
cells. 

In this paper, we present a form of reflective confocal interference microscopy [25] based 
upon laser-feedback interferometry [26-30] combined with phase-shifting techniques which 
yields the phase and amplitude of the reflected field [29-31]. After presenting the basis of the 
method, we apply the technique to measure integrin adhesions in both a fixed and a motile 
cell. Although the morphology of these adhesions has been characterized, the mechanisms that 
govern the formation and lifetime of these adhesions remains a hot topic of research [32,33]. 

In order to relate the phase and amplitude of the reflected field that we measure to the topog- 
raphy and reflectivity of cellular features near the ventral plasma membrane, we have developed 
an algorithm that can be used to separate the contribution of the elastically scattered light by 
sub-cellular structures (or a calibration sphere) from the reflection at the coverslip-buffer inter- 
face. In essence, we solve this inverse problem [34] so that we can find the local topography and 
reflectivity of the sample from a set of interference measurements. We also address the forward 
problem such that, given a priori information about the shape and index of refraction (e.g. of a 
calibration sphere), we can predict the experimentally observed phase and reflectivity [35, 36]. 

Although we are predominately interested in the solution to the inverse problem because the 
majority of applications lack a priori information regarding the object's reflectivity and shape, 
the solution to the forward problem is useful for interpreting data for calibration. As an appli- 
cation of the forward problem, the results of two closely related experiments are analyzed: (1) 
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linear translation of a planar coverslip-buffer interface away from the focal plane of a micro- 
scope objective and (2) reflection from the surface of a polystyrene sphere < 100 nm above 
the coverslip-buffer interface. These two experiments provide a calibration of the profile of the 
focused light from the objective and they lay the foundation for a comparison with results of 
our inverse method. 

Subsequently, we demonstrate that the inverse method as applied to phase data collected 
from the sphere is consistent with the solution of the forward problem. Finally, we apply our 
inverse method to phase data acquired from light that was reflected from stress fibers and focal 
adhesions on the ventral plasma membrane, which lies in apposition to the coverslip-buffer 
interface. In contrast to the application of RICM, where a high concentration of protein in focal 
adhesions may appear black [37], phase-shifted laser feedback interferometry combined with 
our solution of the relevant inverse problem clearly demonstrate that these regions provide a 
strong, bright signal indicating that they have higher reflectivity than the surrounding areas. 

2. Methods 

2.1. Instrument 

As illustrated in the schematic in Fig. 1, the scanning laser feedback interference microscope 
consists of a custom-designed interferometer coupled to an inverted microscope (TE2000U, 
Nikon). The interferometer is based upon a low power, continuous-wave (CW) helium-neon 
laser (1107P, Uniphase) whose output light passes through a linear polarizer (5524, New Focus) 
and a broadband electro-optic phase-modulator (4002, New Focus); the polarization is aligned 
such that only pure phase modulation is achieved. After passage through a subsequent linear 
polarizer, the laser beam is expanded with a beam expander so that the TEMoo mode slightly 
under-fills the back aperture of a microscope objective (CFI60 lOOx, NA =1.45 or a 60x, NA 
=1.3, Nikon). The light enters the microscope via the left-side port and illuminates the sample 
without passing through a tube-lens, a dichroic mirror or a filter. The modulated intensity of the 
laser is monitored with a photodetector (1201, New Focus) that collects the small percentage 
of light that is transmitted through the back laser mirror; the voltage from the detector was 
sampled using an analog-to-digital convertor (DAS 1802, Keithley). 

In order to implement the phase shifting algorithms, five discrete voltages were sequentially 
superimposed upon the modulator using an analog-to-digital convertor (DAS 1802, Keithley) 
and high-voltage function generator (3211, New Focus). In order to build up an image from 
a series of discrete measurements, the sample was translated using a closed-loop piezoelectric 
stage that can move independently along three orthogonal axes (LP100, Mad City Labs); x, y 
are in the transverse plane and z is along the optical or axial axis. Both the spacing between 
pixels (pixel step size) and the direction of the scan were variable. For the axial scans (Section 
4.1), the step-size along the optical ("z" axis) was set to 20 nm; ten scans were acquired and 
averaged. When imaging either an interface or spheres, the step-size along either of the two 
transverse axes was set to 100 nm between pixels. When imaging cells, we used a small step- 
size of 50 nm between pixels. The pixel dwell time, which includes the time to acquire the 
raw data with on-the-fly presentation of the raw phase and visibility to the computer monitor 
is w 0.01 seconds; this does not include the post-processing time to implement the algorithms. 
Therefore, the 400 pixel x 600 pixel image covering 20^m x 30/im in Figs. 6(C) and 6(D) with 
50 nm/pixel, required approximately 40 minutes to acquire. 

All samples (in 35 mm microwell glass-bottom dishes (MatTek Corporation)) were held in 
a custom designed tip-tilt insert that fit into the piezoelectric stage. In order to correlate the 
pixels in the interference image with those pixels in the fluorescence image collected with a 
CCD (Coolsnap, Princeton Instruments), a single affine transformation (consisting of a scaling 
and a rotation) was applied to the scan region and the edges of the transformed images were 
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subsequently cropped by a few pixels to produce rectangular regions. Custom written soft- 
ware (Labview, National Instruments) was used to control the equipment. Data analysis was 
performed with either Matlab (Mathworks) or Mathematica (Wolfram Research). 




microscope objective 
microscope body 
beam expander 




computer 



Fig. 1 . Schematic diagram of experimental configuration of the laser feedback interference 
microscope. Linearly polarized light from a low power continuous-wave helium-neon laser 
passes through a broadband electro-optic phase-modulator and is subsequently expanded 
so that the TEMnn mode fills the back aperture of a high numerical microscope objective. 
The modulated intensity due to laser feedback of light from the sample is monitored with 
a photodetector at the back mirror. A single computer controls the phase shifts to the mod- 
ulator, reads the photodetector signal and controls the piezoelectric stage that moves the 
sample. 

2.2. Assays 

We imaged the coverslip-buffer interface above an unmodified #1.5 coverglass of 35 mm mi- 
crowell glass bottom dishes (MatTek Corporation) filled with a-MEM media (Invitrogen). We 
used immersion oil (Cargille type DF; Cargille Laboratories) between the objective and the 
glass coverslip. Based upon data supplied by the manufacturer, we calculated the index of re- 
fraction at the imaging wavelength (A = 0.6328/im) to be n = 1 .515. We performed calibration 
experiments using polystyrene spheres that had 15/im nominal diameter (FocalCheck, Molec- 
ular Probes). These spheres were placed in a 35 mm microwell plate that contained water and a 
Poly-L-lysine coated (P8920, Sigma, diluted to 0.01%) coverslip. The majority of spheres set- 
tled at a small unknown distance, h 0 above the coverslip such that Brownian motion was essen- 
tially undetectable. From the Cauchy dispersion relation, we determined the index of refraction 
for the polystyrene sphere at the imaging wavelength (A = 0.6328/im) to be n = 1.582 [38]. 



#146921 - $15.00 USD Received 2 May 2011; revised 5 Jul 2011; accepted 15 Jul 201 1; published 27 Jul 201 1 
(C) 201 1 OSA 1 August 201 1 / Vol. 2, No. 8 / BIOMEDICAL OPTICS EXPRESS 2422 



2.3. Cell Culture 

We imaged both fixed and live non-metastatic rat mammary adenocarcinoma MTC cells. Cells 
were cultured in a-MEM media supplemented with 5% FBS and pen/strep (Invitrogen) [39,40]. 
Cells were trypsinized and plated onto the microwell dishes 24 hours prior to fixation and 
staining. The cells were starved for one hour using bovine serum albumin and then stimulated 
by epidermal growth factor for one minute. Finally, cells were rinsed two times with PBS and 
fixed with 4% paraformaldehyde diluted in PBS before being permeabilized with 0.1% triton 
X-100/PBS for 10 minutes. Focal adhesions were identified using indirect immunofluorescence 
with a primary antibody against paxillin (clone 349) (BD Biosciences) and a goat-anti-mouse 
IgG secondary antibody conjugated with Alexa Fluor 488 (Invitrogen). F-actin was stained 
using rhodamine-phalloidin. These probes were easily excited and spectrally separated using 
readily available filter cubes (e.g. FITC and TRITC). 

2.4. Phase Shifting Interferometry 

The basis of phase shifting interferometry is to use an over-determined set of independent 
intensity measurements, so that it is possible to solve for V and O and to eliminate the variable 
I 0 inEq. (1): 

Ii[x) = I 0 {x) {1 +V(x) cos (O(x) + y/,)} (1) 

We employed a popular algorithm that uses five discrete phase shifts in order to reduce the 
sensitivity to random errors in the phase shift, [41]: Xj/j = { — 7t, — 1,0, f ,^}- The visibility and 
phase can then be determined from the measured intensities using Eqs. (2) and (3): 

2{I x +h+2h+h+h) 



and 



2{I 2 -h) 



0 = tan-'( —, — -^-) (3) 

The phase, O, in Eq. (3) is determined from the inverse tangent function such that — § < 
O < j. In order to map the phase to a range 0 < O < 2n, the sign of the numerator and 
denominator in Eq. (3) are examined separately and a factor of n is either added or subtracted 
to the phase [42,43]; the closely related, two argument function "atan2" defined between — it < 
<&<n automatically accounts for the quadrant for the phase and accomplishes the same result. 
Because the phase is determined modulo 2n, visible 2n discontinuities between pixels may 
occur which may be removed using unwrapping algorithms that, depending upon the sign of 
the phase discontinuity between adjacent pixels, either add or subtract an integer multiple of 
2k. 



3. Modeling the Visibility and Optical Path Length 

As indicated in Section 2.1, the basis of our phase-shifted laser feedback interference (LFI) 
imaging is an inverted microscope where sample illumination is provided by a low-power, 
continuous-wave, linearly polarized laser. Using a piezoelectric stage that translates the sample, 
we produce an image by collecting data line-by-line in a two-dimensional raster format with 
a scan line that begins in a "sample-free" region that contains only buffer above the coverslip 
and proceeds pixel-by-pixel into a "sample" region that contains either cells or a calibration 
sphere. At the end of the scan-line, the scan position is returned to the "sample-free" region and 
is incremented along the orthogonal direction. Unless we specifically want to characterize the 
effect of defocus, the coverslip-buffer interface is held fixed at the focus of the objective and an 
axial scan is not performed. 
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Using phase-shifted LFI microscopy (psLFIM), the steady-state change in the intensity of 
the laser caused by weak feedback of elastically scattered light may be written in a form 
that is commonly used to represent the interference of two coherent waves [27,29]: = 
I 0 (x) {1 + V(x) cos (O(x) + y/j)} where I u (x), V(x) and <S>(x) represent: a background inten- 
sity; the Michelson fringe visibility and the optical phase difference between the two fields, 
respectively. After introducing a discrete, experimentally controlled phase shift, y/j, the corre- 
sponding intensity measurement, /,, is measured (henceforth, we will suppress the notation that 
represents the spatial variation of the parameters). The phase shifting algorithm can be applied 
to the measured intensities, in order to determine the visibility and phase independently of 
one another and independently of the background intensity, I 0 . 

In LFI, the fringe visibility, V, is proportional to the electric field that is reflected by the sam- 
ple; because the coupling efficiency of the laser feedback interferometer is less than 100%, the 
fringe visibility is always less than the maximal value that would be calculated purely on the ba- 
sis of the Fresnel reflection coefficients (previous work indicates a coupling efficiency of about 
40% [27-29]). Because the fringe visibility is intimately related to the reflection coefficient, we 
will now use the terms "fringe visibility" and "reflectivity" interchangeably. 

As with other forms of interferometry, the change in optical phase is related to the optical 
path length, OPL, between the points pi and p 2 , as: O = x 0PL = T Ipi ne ^ s f 44 J- Near the 
focal plane of a high numerical aperture objective, however, the phase changes nonlinearly 
with defocus [45,47,48] and accumulates a phase shift of — 5 for positive defocus. When 
the sample is imaged close to the focal plane, the effect of the nonuniform spacing of the 
interference fringes may be approximated with a constant factor, / (> 1), so that the optical 
path length is related to the phase as: OPL = hf® [45,47]. 

Regardless of the location of the scan position, we can interpret the measurement as con- 
taining a contribution from "signal" and "background". In the "sample-free" region, the re- 
flected light that re-enters the laser consists of a reflection from the coverslip-buffer interface 
("signal term") and a considerably weaker reflection from the optical elements along the path 
("background term"); because of the index-matching oil, the reflection from the bottom of the 
coverslip is negligible. These two reflected fields give rise to an observable interference signal. 
After translating the scan position to the "sample" region, we interpret the "signal term" as 
consisting of light reflected from the sample (cell or calibration sphere) which interferes with 
the "background term" that arises from the field reflected from the coverslip-buffer interface 
and the reflections from the optical elements. 

In either the "sample-free" or "sample" region, we model the interference from the multiple 
reflections as arising from two discrete fields that interfere inside the laser cavity. These multi- 
ple reflections may be modeled as a change in the effective reflectivity of the laser mirror such 
that the measured intensity may be written as Eq. (4) [26,28,29]: 

h = I 0 {1 +Acos (0 + \f/j) +Bcos (</> + 6 + yi)} ( 4 ) 

In the "sample-free" region, the parameters {A, 0} in Eq. (4) represent the visibility and ac- 
cumulated phase associated with the weak reflections from the optics along the path ("back- 
ground term") and the parameters {B, 0} represent the subsequent reflection that occurs after 
the partially transmitted field is reflected from the coverslip-buffer interface ("signal term"). By 
contrast, when the scan position has been translated into the "sample" region, the visibility A 
represents a "background term" which consists of the reflection from the coverslip-buffer inter- 
face (and the weaker reflection from the optical elements) and B represents the "signal term" 
that arises from light that is reflected from the sample. Accordingly, in either region, (j) and 6 
represent the phases that have accumulated following the reflection from the primary interface 
(optical elements or coverslip-buffer) and the secondary interface (coverslip-buffer or sample), 
respectively. 
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A little algebra allows Eq. (4) to be represented as Eq. (5) which is then amenable to the 
application of the phase-shifting algorithm: 



4=/ 0 {l+mcos(y+yi)} (5) 
where the visibility, to, and phase, y, are related to the variables in Eq. (4) as: 



to = yA 2 +fi 2 +2ABcos(0) (6) 

and 

_ u fisin(0) . , 

In the "sample-free" region, determining the parameters {B, 6} associated with the reflec- 
tion from the coverslip-buffer interface is primarily relevant. Once the scan position has been 
translated to the "sample" region, however, we wish to determine the parameters {B, 6} which 
are associated with the reflection from the sample. 

3.1. Predicting the Visibility and Optical Path Length Using A Priori Information: the For- 
ward Problem 

We formulate the forward problem: given estimates of the fringe visibility and phase {A, (j)} and 
{B, 6}, predict the measured values m and y. Using Eqs. (6) and (7) and a priori information 
so as to estimate {A, 0} and {B, 0}, we make an initial prediction of the visibility and phase, 
m and y, and refine the estimate after comparison with the measured values. 

3.2. Determining the Topography and Index of Refraction of a Sample Without A Priori In- 
formation: the Inverse Problem 

As an alternative to the forward problem, it is possible to formulate the inverse problem: given 
measurements of the visibility and phase {m, y}, find the unknowns {B, 0} associated with the 
reflection from the sample. The approach to the inverse problem as applied to reflections from 
multiple interfaces begins with a direct solution of Eqs. (4) and (5) for the reflectivity, B, and 
phase, 6: 



B= ^m 2 +A 2 -2mAcos(y-<$>) (8) 

and 

e=*n-H m 7 (y :f J - (9) 

mcos(y— q>) —A 

Therefore, provided that the values {A, 0} have been measured, Eqs. (8) and (9) may be 
combined with a measurement of {to, y} in order to determine the value of the visibility and 
phase, {B, 6}, associated with the reflection from the sample. 

4. Results 

4.1. Application of the Forward Problem: Change in the Visibility and Optical Path Length 
Following Defocus of the Coverslip-Buffer Interface 

Figure 2(a) and 2(b) (dotted lines) shows the visibility, to, and measured position of the inter- 
face, Az m (obtained from Az m = = 4^7 y with X = 0.6328/im and n = 1.515) measured in 
the "sample-free" region as a function of the imposed displacement, Az, of the coverslip-buffer 
interface. The data in Fig. 2(a) show that as the defocus increases, the visibility oscillates with a 
superposed decay envelope. The measured stage position, Fig. 2(b), appears to increase linearly 
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and there are superposed weak oscillations that become increasingly pronounced with defocus. 
Although a departure from linear phase variation with a small amount of defocus is not readily 
evident in Fig. 2(b), this effect is clearly observed after a line with slope Ann/X is subtracted 
from the data. Fig. 2(c) (dotted lines) shows Ay = y— (Ann/X) Az as the coverslip-buffer inter- 
face was translated away from focus by Az = 2jj.m. 

In order to predict the observed visibility and phase data using Eqs. (6) and (7), we need 
initial estimates for A, B and 6; without loss of generality, the overall phase factor, 0, may 
set to zero. Initially, we estimate the maximum reflection from the interface as B = fc|r|, where 
K = 44% is the expected coupling efficiency and r is the Fresnel coefficient at normal incidence; 
using the index of refraction for glass and buffer as 1.515 and 1.333, respectively, we have 
B rj 0.44 x 0.064 s» 0.028. It is more challenging to determine an initial estimate for A, however, 
if we assume that 2% of the incident field is reflected from the optical elements (i.e. a reflection 
that is consistent with the anti-reflecting coating at the back of the objective), then we can 
estimate A « 0.44 x 0.02 x 0.20 s» 0.0018. 

We model the observed envelope of the visibility in Fig. 2(a) by assuming that the reflection 
from the coverslip-buffer interface falls off with defocus as the square-root of a Lorentzian 
function; this functional form is a good approximation to the axial profile of the electric field 
of our Gaussian-Lorentzian TEMoo laser beam focused by the microscope objective [45,46]. 

Therefore, we have B = B ( ,j \j 1 + (^) 2 where y/3o is the half-width at half-maximum, B a = 
0.028, G sw 0.3jUm for our microscope objective [45] and Az is measured in microns. 

The linear phase change in Fig. 2(b) is proportional to the imposed linear translation, Az, of 
the coverslip-buffer interface as 9 = (2n/X) «2Az where n = 1.515 and A = 0.6328jUm. There- 
fore, the analysis of the forward problem may be used to predict the measured visibility and the 
change in phase using Eqs. (6) and (7), respectively. Figures 2(a) and 2(b) (solid lines) show the 
result of a non-linear least squares fit with parameters A, B and G (with initial estimates of A = 

0.0018, B = 0.028 and a = 0.3, respectively) such that A = 0.0019 andfi = 0.028/yi + (g%) 2 

with 0 = 30.00 rad^m-'Az. 

Theory predicts that the observed phase also depends upon the numerical aperture of the 
microscope objective [45]. Because of the nonlinear variation with defocus (Fig. 2(c)) there 
is an accumulated phase difference —n/2 at large defocus which is known as the Gouy 
phase shift [48]. Figure 2(c) (solid black line) shows this predicted phase anomaly, Ajo = 
— tan -1 (2 Az/zo), where z ( , = 0.3 /im. Therefore, it is appropriate to refine our predicted linear 
phase variation, 6 = (Ann/X)Az, by including a nonlinear term — tan -1 (2Az/z 0 ). Figure 2(c) 
(solid red line) shows Ay based upon Eq. (7) using 6 = (4nn/X)Az — tan -1 (2 Az/z„) + n 
and the parameters that resulted from the fit to the data in Fig. 2(b) (i.e. A = 0.0019; B = 
0.028 /Vl + (Az/0.28) 2 ; n = 1.515 and X = 0.6328 /im). Following convention, we express 
the variation of phase near the focal plane for small values of defocus as 9 « (4nn/X)f~ 1 Az 
where / _1 is a constant factor [47]. Using a linear expansion of tan -1 (2 Az/z 0 ), we obtain 
f~ l i=s 0.8. Conversely, given a phase measurement near focus, 0, the interface is at height 
h rj above the focal plane, where / = 1.2. 

4.2. Application of the Forward Problem: Reconstruction of the Visibility and Optical Path 
Length to the Surface of a Sphere at a Small Distance from the Coverslip-Buffer Interface 

In a second closely related demonstration of the application of the forward problem, we imaged 
polystyrene spheres and measured the interference between light reflected from the sphere and 
light reflected from the coverslip-buffer interface. This relatively large sphere was imaged with 
the focus of the objective held fixed at the coverslip-buffer interface; this approach assured that 
the reflection from the (bottom) surface of the sphere closest to the coverslip-buffer interface 
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Fig. 2. Effect of defocus on visibility of the interference fringes from a coverslip-buffer 
interface, (a) Change in the measured visibility, m, (dotted lines) as the coverslip-buffer 
interface was translated by Az; solid line is based upon a solution to the forward problem 
(Eq. (6)). (b) Position of the interface, Az m = X y/(4jrn), as a function of defocus (dotted); 
solid line is based upon Eq. (7). (c) Difference between measured and predicted phase, 
Ay = y — (4m/X.)Az, near focus (dotted). Plot of Aye = — tan -1 (2Az/z 0 ) (solid black 
line) where z 0 = 0.3 (lm. Predicted variation (solid red line) on the basis of Eq. (7) with 
6 = (4m/ 1) Az - tan" 1 (2Az/ Zo ) + n. 



was greater than the (top) reflection from the surface of the sphere that was furthest from the 
interface. 

By analogy with the analysis of the planar interface, we address the forward problem: given 
the fringe visibility and phase associated with the reflection from the coverslip-buffer interface 
{A, 0} and the sphere {B, 6}, find m and y on the basis of Eqs. (6) and (7). From the initial 
experiment with a planar interface, we are able to estimate the reflectivity, A, from the value 
of the fringe visibility when the interface is in focus, i.e. Az = 0 in Fig. 4(a), as A = 0.035. In 
addition to representing the dominant reflection from the coverslip-buffer interface, the term A 
also accounts for the constant, considerably weaker reflections from the optical elements. 

When imaging the reflection from a tilt-free, planar reflector at the focus of a high numer- 
ical microscope objective, all of the incident rays will be collected by the objective, however, 
for a sphere the angle of reflection of an incident ray at the sphere's surface depends upon the 
transverse scan position. Therefore, unlike the reflection from the planar interface, the angle 
of reflection of an incident ray at the sphere's surface depends upon the transverse scan po- 
sition and the visibility, B, and OPL depend upon the angle of incidence. Because we do not 
have an analytical description for the bidirectional reflectance distribution function (BRDF) 
that describes the reflection from the sphere [49], we approach the fit to the data by invoking a 
"tangent-plane approximation" which assumes that light is reflected from a diffraction-limited 
patch on a local tangent plane at the sphere's surface [50]. Therefore, when the scan position is 
directly beneath the center of the sphere, light is reflected from a tangent plane that is parallel 
to the plane of the coverslip-buffer interface, and consequently the optical path will be parallel 
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to the optical (z) axis. 

As shown in Fig. 3, when the scan position is not at the center of the sphere, the majority of 
incident rays that reflect from the sphere's surface (blue lines) are reflected at an angle such that 
they do not re-enter the objective (red lines). Instead, only those rays that are reflected from a 
diffraction-limited patch on a local tangent plane at the sphere's surface (shown in black) are 
collected by the objective. Therefore, as the scan position moves outward from the center of 
the sphere, the optical path can be determined by tracing the path of the ray between the focus 
and the point where the tangent plane intersects the sphere [51]. Therefore the retro-reflected 
ray follows a path which forms an angle, with respect to the optical axis. 




Fig. 3. Imaging geometry when scanning a sphere. When the scan position is directly be- 
neath the center of the sphere, light is reflected from a tangent plane that is parallel to the 
plane of the coverslip-buffer interface, consequently the optical path will be parallel to the 
optical (z) axis. As the scan position moves outward from the center of the sphere, the op- 
tical path can be determined by tracing the path of the ray from the microscope objective 
that intersects the sphere normal to the tangent plane and reflects back into the objective. 
Therefore, the optical path, 8 = 0PL/(2n), at the scan position, x and the perpendicular 
distance, h, from the interface (at x 7 ) to the tangent plane. 



We define the center of the scan as the unique position where light is reflected from a tangent 
plane that contains the coverslip-buffer interface such that the optical path will be parallel to 
the optical (z) axis (Fig. 3). At this location, OPL = n2h 0 where h 0 is the height of the sphere 
above the coverslip-buffer interface and n is the index of refraction of the buffer. Therefore, at 
the transverse scan position (x,0), the retro-reflected ray from a sphere of radius R follows a 
path which forms an angle, #, with respect to the optical axis given by Eq. (10): 

t> = tan- 1 (-^-). (10) 
R + h 0 

The retro-reflected ray that follows a path that makes an angle # with respect to the optical 
axis, travels a distance 25 given by Eq. (11): 



5 = ^(R + h 0 ) 2 +x 2 -R. (11) 

Correspondingly, the OPL = n2 8 and the phase is given by: 

2k 

e = -r-0PL + 7r. (12) 

A 

Figure 4(A) shows the measured fringe visibility, m, obtained using the phase shifting algo- 
rithm (Eqs. (2) and (3)) applied at a series of discrete scan positions covering a square; the scans 
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positions were spaced at IOOj^j and we used a 60x, NA =1.3 (Nikon) objective for imaging. 
Figure 4(B) (dotted line) shows the visibility along a single row of pixels corresponding to 
y = 0 and —6/im < x < 6/im. Correspondingly, the measured phase, 7, (via Eq. (3)) is shown 
in Fig. 4(C) (dotted line). Although the microscope focus is held fixed on the coverslip-buffer 
interface, the visibility and phase data show oscillations with a superposed envelope that decays 
as the scan position moves outward from the center of the sphere. 

In order to fit the experimental data, we account for the variation in the density of rays 
with scan position because at the center of the scan, all of the incident rays contribute to the 
interference, but as the scan location moves off-center, fewer of the paraxial rays reflect from 
the sphere and re-enter the microscope objective. We introduce an angular dependent term and 
express the reflectivity from the sphere's surface as: B = Boe^^l®"" 1 ' / \J\ + (h/a) 2 , where B„ 
is the fringe visibility at the center of the scan (i.e. # = 0) and & 0 is a fitting parameter. 

Using an initial estimate B a = 0.038, based upon the index of refraction of the sphere and 
buffer as n s = 1.582 and n = 1.333, respectively and a coupling efficiency, K = 44%, we fit 
the data in Figs. 4(A) and 4(B) with Eqs. (6) and (7). The fits to these data are shown in Figs. 
4(A) and 4(B) (solid lines) using: A = 0.0345; B = 0.045 e'^^ 2 /y/l + (h/a) 2 ; A = 25"; 
a = 0.3/im; R = 7 '.7//m; h () = 0.005/im; n = 1.333 and A = 0.633/im; the values for the index 
of refraction, n, and the wavelength, A, were held fixed. Consequently, we can now use the 
values A = 0.0345 and B 0 = 0.045 to refine the initial estimate of the index of refraction of the 
sphere. If we assume that the reflected field from the coverslip-buffer interface couples back 
into the laser with the same efficiency, K, as the field that is reflected from the buffer-sphere 
interface, we can eliminate the dependence upon K and we can use the ratio B 0 /A to determine 

"sphere = 1-58. 

4.3. Application of the Inverse Problem: Reconstruction of the Shape of a 15/im Diameter 
Sphere Using the Inverse Method 

We now demonstrate the feasibility of approaching the inverse problem such that we are able 
to recover the shape and height of the polystyrene sphere held at a fixed distance above the 
coverslip-buffer interface without a priori information that we are imaging a sphere. Specifi- 
cally, we wish to find the visibility and phase {B, 0} and reconstruct the shape and height of 
the surface above the interface from measurements of the visibility and phase {m, 7}. 

Equations (8) and (9) relate the visibility and phase {B, 6} associated with the reflection 
from the unknown object to the measured visibility and phase {m, 7} and {A, 0}associated 
with the reflection from the coverslip-buffer interface. These two parameters may be determined 
by collecting visibility and phase data from the coverslip-buffer interface in the "sample-free" 
region where B equals zero. In this region, Eqs. (6) and (7) combined with a measurement of 
{m, 7} yield {A, (j)}. With A and 0 determined in the "sample-free" regions at the beginning 
and end of each scan-line, the values of {A, (j)} in the "sample" region may be determined by 
interpolation. 

Figure 5(A) shows the visibility, B, obtained by applying the inverse method to the raw 
visibility, m, and phase data, 7. Unlike the data in Fig. 4(A), however, reconstruction of the 
visibility using the inverse method eliminates the series of interference fringes. Because the 
inverse method reconstructs the visibility and phase associated with the reflection from the 
sample, the "sample-free" regions are shown as black. From Fig. 4(A), the measured reflectivity 
is highest in the center of the bead and lower towards the periphery of the scan region; this 
result is a consequence of the combined effect of the axial profile of the incident illumination 
and the curvature of the surface of the sphere. The visibility data taken from a single row of 
pixels corresponding to y = 0 and —6/im <x< 6/im in Fig. 5(A) (which corresponds to the 
same pixels as Fig. 4(B)) are shown in Fig. 5(B) (dotted lines). The measured OPL, calculated 
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Fig. 4. (A) Visibility of the fringes from the surface of a 7.7 jj.m radius polystyrene sphere 
in water. (B) Fringe visibility and (C) phase from a single row of pixels; the analytical fits 
(solid lines) are based upon analysis of the forward problem and dotted lines are experi- 
ment. 



from Eq. (12), as a function of scan position, x, is shown in Fig. 5(C) (black squares). From 
the data, the height of the center of the sphere above the interface, h 0 , was determined to be 
h a = 0.005/im. The 0PL data were compared to the predicted value using Eq. (11) (solid black 
line) and R = 7.7jUm. 

After the application of this inverse method, the phase data, 6, showed discrete phase jumps 
and conventional phase unwrapping algorithms were sufficient to unwrap the phase yielding the 
data shown in Fig. 5(B). Therefore, it was possible to reconstruct the 0PL from the unwrapped 
phase, 9 as: 0PL = ^[d-7t\. 

In order to verify that the inverse method reproduces the expected relationship between the 
height h and 0PL, we compared the measured height of the sphere above the coverslip-medium 
interface, h, as a function of the distance from the center of the sphere, x 1 with the expected 
value (i.e. h = h(x')). From the geometry of Fig. 3, at the scan position, x, the shortest distance 
between the local tangent plane to the sphere's surface and the coverslip-buffer interface, h, is 
given by Eq. (13) which is measured at a virtual location, xf, given by Eq. (14): 

h = 5cos(?>), (13) 
x = x-5sm(-&), (14) 

where # and 8 are given by Eqs. (10) and (11), respectively. 

Figure 5(C) (red squares) shows the height of the sphere above the coverslip-medium inter- 
face, h, as a function of the distance from the center of the sphere, x! determined using Eqs. 
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Fig. 5. Inverse method used to determine the visibility and height of the surface of a sphere 
above an interface. (A & B) visibility and (C) 8 (black squares) and height (red squares) as 
a function of the virtual scan position, xf (h = h(x')). The fit to 8 (black line) and predicted 
height of sphere, h = h(x') (red line) is also shown based upon a sphere with radius R = 
7.7/xm and h„ = 0.005/im above the interface. Fit to visibility using Eq. (8) (solid line in 
B). 



(10)— (14). The re-mapped data were then compared to the predicted value (solid red line) for 
a sphere with radius, R = 7.7/im. Using h = h(x'), we compared the visibility data, m, to a 
prediction using Eq. (8) with B = 0.045 e^^fy/l + (h/o) 2 with = 25°; G = 0.3/im 
(Fig. 5(B), solid line). 

5. Imaging the Reflectivity and Topography of the Ventral Surface of the Cell 

5.1. Reconstruction of the Visibility and Optical Path Length of Stress Fibers and Focal Ad- 
hesions 

We now apply the inverse method to phase data that was acquired from a fixed MTC cell plated 
on glass bottom dishes. This approach reconstructs the visibility and phase associated with re- 
flective sites near the ventral plasma membrane which lies in apposition to the coverslip-buffer 
interface. In order to gain insight into the location of dense, clustered sub-cellular structures 
and their relationship to the interference images, we used indirect immunofluorescence to im- 
age both F-actin and paxillin in a single cell (Figs. 6(A) and 6(B), respectively); the 20 /im x 
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30 jUm scan region also shows black portions of the coverslip to the left of the periphery of the 
cell. The bright, dense sites of the adhesion-associated protein paxillin in Fig. 6(B) indicate the 
presence of focal adhesions [52] and the F-actin distribution in Fig. 6(A) may be interpreted as 
either ventral and/or dorsal stress fibers; these dense bundles of actin filaments attach to focal 
adhesions with the ventral fibers aligned nearly parallel to the ventral surface and the dorsal 
fibers aligned toward the dorsal surface with one of their ends tethered to the adhesions [53]. 




0 2 4 6 8 10 12 14 16 x1 ° 3 0.1 0.15 0.2 0.25 0.3 0.35 \im 



Fig. 6. Fluorescence and interference images obtained near the ventral plasma membrane of 
fixed cells. (A) Immunofluorescence image showing F-actin localization near the periphery 
of the cell, (scale bar: 5jJm). (B) Immunofluorescence image of the density of paxillin. 
(C) Fringe visibility and (D) height of the reflective features above the coverslip-buffer 
interface at discrete scan points that cover the same region as the fluorescence images 
calculated using the inverse method (Eqs. (8) and (9)). The height, z, was determined from 
the phase, 6, as: z = 9 with A = 0.6328(tm, n h = 1.333 and / = 1.2. Using this 

method, the height map represents the distance above the coverslip. 

The application of the inverse method (Eqs. (8) and (9)) to data acquired in this 20 /im x 
30 jiim region yields a map of the visibility (Fig. 6(C)) and the height (Fig. 6(D)), respectively. 
The color map in Fig. 6(C) shows that the visibilities range from 0.005 to less than 0.02 and 
that the majority of scan points have very low visibility. From the data in Fig. 6(C), it may 
be observed that in the "sample-free" region to the left of the cell, the application of Eq. (8) 
produces B essentially equal to zero. Ideally, Eq. (8) would yield B = 0 when the scan region 
is on the "sample-free" region because the measured visibility, m, and phase, y, are: m = A and 
y = (j). It may be observed, however, that for some pixels, B ss 0. This deviation results from 
using a linear fit to the measured values of A and 0 along an entire row of pixels rather than 
the precisely measured values of A and 0 . Therefore, the small departure of B from zero in the 
"sample-free" region indicates that there is a small error resulting from the linear fit. 

By comparison, under the cell, there are clusters of pixels with relatively high visibility and 
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the location of these bright regions are observed to be correlated with either the regions of high 
F-actin density shown in Fig. 6(A) or the paxillin-dense regions shown in Fig. 6(B) or with both 
F-actin and paxillin. 

From the visibility, it is possible to estimate the Fresnel coefficient at each scan location 
(assuming that the angular variations of these coefficients are negligible) and consequently the 
index of refraction at each pixel. Therefore, Fig. 6(C) may be interpreted as a map of the spatial 
variation of the index of refraction with scan position. The visibility, B, obtained from Eq. 
(8) is related to the coupling efficiency, K", and the Fresnel coefficient at normal incidence, r, 
as B = fc|r|. In the region "sample-free" region outside the cell, Eq. (8) predicts that B = 0 
at the coverslip-buffer interface. When the scan position is located beneath a region of the cell 
that is predominately cytoplasm (with index of refraction equal to 1 .35), we expect the visibility 
B rs 0.44| 35 I = 0.003. From the data, the majority of the pixels have an index of refraction 
1 .33 < n < 1.35, but the higher visibility pixels (those with yellow to red values on the visibility 
colormap), have an index of refraction such that 1 .4 < n < 1 .43. 

The height map in Fig. 6(D) is determined, as explained in Section 4. 1 , by assuming that light 
propagates in a straight line from the objective through the buffer (i.e. see the geometry of Fig. 
S4 when the scan position is at the center of the sphere) until it was reflected from a discrete 
site near the ventral membrane. Therefore, the height, z, was determined from the phase, 0, 
as: z = ffj^ 9 with A = 0.6328/im, ny = 1.333 and /= 1.2. In the "sample-free" region, the 
value of 6 is essentially meaningless because the visibility B should be identically zero in this 
region. Nevertheless, the application of Eq. (9) yields a value at each pixel in the field-of-view 
including the pixels in the "sample free" region. It may be observed from Fig. 6(D), that at the 
pixels close to the cell edge, the heights are essentially constant, but that some rows of data 
have a slight variation as the scan position moves to the edge of the "sample free" region. This 
variation results from the use of a linear fit to the measured values of A and (j) in this region 
rather than the precisely measured values. 

6. Visibility and Topography at the Trailing Edge of a Motile Cell 

Figure 7 shows six scans of the visibility and topography (top row and bottom row, respectively) 
near the edge of a live cell. The data for each scan were acquired in approximately 25 seconds 
and the total elapsed time to acquire the six scans was approximately 40 minutes. As in Section 
5. 1, the data were reconstructed using the inverse method such that the visibility and topography 
represent the actual visibility of the cellular features and the topography represents the height 
of the reflective structures above the glass. Unlike the data obtained from fixed cells, however, 
neither paxillin nor F-actin were labeled in this cell. 

In each of the six data sets, the dark blue regions (essentially zero visibility) represent the 
glass substrate. Examination of either the visibility data (top panel) or topography, indicate that 
the bulk motion of the cell is toward the south-east (bottom-right) corner. The live cell data 
clearly demonstrate dynamic remodeling of the cytoskeleton and the dynamics of the trailing 
edge. In each of the first five scans, three filopodia may be observed. By the sixth scan, essen- 
tially all of the cell has disappeared from the field-of-view. 

In the first scan, two filopodia appear to have approximately uniform diameter (« 0.5/im) 
and a third considerably thicker and nonuniform filopodium may be observed. Between the 
three filopodia of the first two scans, the edge of the cell is visible. This region has a visibility 
w 0.006. By comparison with the fluorescence and visibility data shown in Fig. 6(A) and 6(C), 
the data from this edge is consistent with the thin dendritic actin network of the lamella. The 
visibility data from the first scan also clearly indicate the presence of highly reflective material 
in the filopodia (with B « 0.01). In two of the three filopodia, this reflective material is most 
densely located at the mid-section of the shaft of the filopodia and at the base (the region 
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where the filopodium meets the lamella), whereas the thick filopodium has very densely packed, 
highly reflective material throughout. Because this material has the same reflectivity as the focal 
adhesions in Fig. 6(C), it suggests that these regions represent the dense aggregation of actin 
and paxillin. The height data in these regions, i=a lOOnm above the extracellular matrix, is also 
consistent with the data from focal adhesions in Fig. 6(D). Therefore, on the basis of both 
the reflectivity and the height data, it may be concluded that these dense regions represent 
adhesions. Compared with the filopodia, the lamella in the first scan appears to be situated w 
160 nm above the matrix which is « 60nm higher than the adjacent filopodia. 

From the data in the first two scans, it may be observed that the transverse speed of the tips 
of the filopodia (ss 66^ ) is much slower than the speed of the retracting lamella (m 250^). 
Although the slow retraction of the tips of the filopodia may be observed throughout the first 
five scans, between the fifth and sixth scan, the filopodia disappear completely from the field- 
of-view. By the sixth scan, only three dense, reflective regions remain; two of these regions 
correlate with the leftmost filopodium and the third blob correlates with the central filopodium. 
Comparing the location of these adhesions with the earlier scans, it may be observed that these 
adhesions were present at the base of the filopodia (a hint of the adhesion from the base of the 
third, thick filopodium may also be observed in the sixth scan). 
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Fig. 7. Visibility and topography reconstructed for a live cell using the inverse method. 
The cell may be observed to be moving towards the south-east (lower-right) of the image. 
Top row: six scans of the visibility. Bottom row: six scans of the topography. Each of 
the scans contains 60 x 40 pixels and was acquired in approximately 25 seconds. The 
temporal separation between scans was about 6 minutes and the total time elapsed time 
was approximately 40 minutes. 



7. Discussion 

We have coupled phase-shifted laser feedback interferometry with an inverted microscope so 
as to measure the change in phase and amplitude of the electric field upon reflection from 
discrete sites near the ventral plasma membrane when the cell was imaged through a glass 
coverslip. Combined with the experimental implementation, we have applied a solution to the 
forward problem appropriate to our imaging geometry so that we were able to compare model 
predictions with the measured phase change near focus and with the reflectivity and phase of 
the field reflected from a calibration sphere. In addition, we have developed an algorithm that 
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addresses the inverse problem so that we can separate the contribution of the reflection from 
the coverslip-buffer interface from the field that has been elastically scattered by sub-cellular 
structures. 

Using this approach, we have reconstructed the topography and reflectivity of regions of the 
ventral membrane at focal adhesions. As expected, there is a strong correlation between the 
presence of paxillin and actin in stress fibers (as identified with fluorescence imaging) and the 
interference images. Unlike the data obtained with RICM, which show dark regions at focal 
adhesions [1^1], using psLFIM we are able to demonstrate that these optically dense regions 
on the cytoplasmic face have high reflectivity and we are able to measure the topography of the 
regions independently of the reflectivity. The data clearly demonstrate that the majority of the 
highly reflective regions are closer to the coverslip than the less reflective regions. 

Although psLFIM offers a high precision measurement of the phase and reflectivity (limited 
by photon counting statistics), an accurate determination of sample topography and index of 
refraction depends upon several model assumptions. In order to interpret our interference data, 
we have assumed that the electric field collected by the microscope objective consists of a re- 
flection of the incident electric field from the coverslip-buffer interface and an additional, single 
reflection from a discrete, protein dense, site near the ventral plasma membrane. Alternatively, 
it is possible to incorporate a model that stipulates reflection from a series of discrete layers of 
known index of refraction [54]. 

Because we have assumed reflection from a single discrete layer, we are able to apply this 
model to the measured fringe visibility and phase and determine the index of refraction and the 
optical path length (topography) at each pixel. From the measured data, we are able to determine 
the distance between the coverslip-buffer interface and the adhesion. This interpretation relies 
on an assumption that actin filaments and the protein plaque are quite dense such that the 
backscattered photons reflect from the surface of the adhesion. 

Assuming a homogeneous index of refraction for the cytoplasm that is essentially equal to the 
index of refraction of the buffer, the reflectivity data indicate that the majority of pixels near the 
ventral plasma membrane have an index of refraction that is close to the index of the cytoplasm, 
but that the index of refraction of sub-cellular sites in focal adhesions have an index greater 
than 1.4. Although this is a precise estimate and the higher index of refraction is consistent 
with measured values from lipid and organelles like mitochondria [55], the accuracy of this 
measurement depends upon knowledge of the index of refraction of the cytoplasm (generally 
assumed to be « 1 .35 [54]) and the validity of the model. 

Furthermore, in order to interpret the reflectivity data at each pixel, we have assumed that 
the incident, linearly polarized light is not significantly depolarized as it undergoes a single 
reflection from discrete sites within the focal adhesions. This assumption is consistent with 
observations of scattering from cells which demonstrated that the dominate sub-cellular features 
that backscatter light, organelles with sizes between A /4 — A /2 (e.g. peroxisomes, lysosomes 
and mitochondria), do not depolarize the linearly polarized incident light [55-57]. Our data 
suggest that this assumption is justified. 

Recently, the fluorescence based, super-resolution technique "interferometric photoactivated 
localization microscopy (iPALM) [58]" has been applied to determine the stratification of pro- 
tein topography within a focal adhesion [59]; this technique offers transverse resolution below 
the diffraction limit [13] and the axial position of the fluorescence protein is determined in- 
terferometrically. Based upon images obtained from fixed human osteosarcoma (U20S) and 
mouse embryonic fibroblast (MEF) cells on fibronectin-coated glass, the data suggest that the 
adapter proteins within the adhesion (e.g. paxillin and vinculin) reside in a thin layer that is 
« 40 nm above the coverslip. By contrast, this method demonstrates that within the focal ad- 
hesion actin is located within a relatively thick layer of width 15 nm at a mean height of 
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pa 96nm above the coverslip. The data from these experiments suggest that the topography that 
we have measured within the focal adhesion, Fig. 6(D), predominately represents the height 
of the actin layer 100 nm) above the coverslip. Therefore, we conclude that the majority of 
the elastically scattered photons that contribute to the most reflective structures in Fig. 6(C) are 
reflected from F-actin rather than from a thin layer of paxillin. 

As shown in Fig. 7, the technique may be used to image live, motile cells. The data demon- 
strate that both the visibility and topography may be determined as the cell rapidly remodels. 
As with other imaging techniques, the speed of the measurements (pixel dwell time) is limited 
by the desired signal-to-noise ratio and by the Nyquist sampling criterion. Because the preci- 
sion of a phase measurement increases as the number of collected photons increases, increasing 
the pixel dwell time will increase the precision. Conversely, rapid motion may lead to blur and 
decrease in signal-to-noise caused by reduce fringe visibility. Previously, we have demonstrated 
that the precision of a phase measurement decreases as the fringe visibility decreases and that 
when the fringe visibility is 0.01, the random error in the height measurement was measured to 
be approximately 1 nm, but the error increases to 10 nm when the visibility equals 0.001 [30]. 

Currently, all of the phase calculations are done on-the-fly and written to screen and plotted. 
There is considerable computational overhead associated with this approach, but it allows the 
user to monitor aspects of the scan. Alternatively, it is possible to implement code that does not 
plot any data as it is acquired. This is approach, coupled with optimal memory management 
of the acquired data, would be faster. Because the total scan time depends upon both the pixel 
dwell and the number of pixels in the scan, there is always a compromise in choosing the pixel 
step size. In order to avoid aliasing, the step size is chosen so as match the effective point spread 
function. Because of experimental noise, however, it is important to sample even faster than the 
Nyquist sampling frequency. The chosen step size therefore depends upon the nature of the rate 
of change of the measured phase, the signal-to-noise ratio and the requisite dwell-time. One 
approach is to implement a coarse scan in order to determine gross changes and then a finer 
scan with smaller step-sizes. 

Applying psLFIM to a live cell can yield data that is extremely difficult to obtain using fluo- 
rescent imaging methods such as total internal fluorescence microscopy. In particular, the data 
obtained from the motile cell in Fig. 7 may be used to gain insight into the physiological mech- 
anisms that govern the retraction of filopodia. Both the reflectivity and topography data (Figs. 
7 top and bottom rows, respectively) suggest that there are adhesions at the tip, shaft and base 
of the two uniform filopodia and this conclusion is consistent with the observations of three 
distinct type of adhesions that have been found in the filopodium of sensory growth cones [60]. 
These three adhesions types have unique functions in filopodium dynamics. In particular, it 
has been shown that: (1) tip adhesions provide guidance cues that govern filopodium move- 
ments; (2) shaft adhesions control the dynamics of the dendritic actin network that constitutes 
the lamella (termed a veil) and (3) basal adhesions affect filopodium dynamics. In the retracting 
filopodia observed in Fig. 7, the tip adhesions appear quite weak whereas the basal adhesions 
appear firmly anchored to the extracellular matrix such that they are ultimately left behind (in 
the sixth scan) as the cell leaves the field-of-view. The extended lifetime of these basal adhe- 
sions is completely consistent with previous observations [60]. The reflectivity and topography 
data associated with the thick filopodium is consistent with previous observations in growth 
cones that suggest that thickened filopodia support strong adhesions, while the uniform and 
even margins associated with the other two filopodia are consistent with previous observations 
associated with moving filopodia. The rapid retraction of the veil region observed between the 
first and second scan indicate that the veil clearly does not adhere to the extracellular matrix. 
This observation is consistent with both the reflectivity and topography data obtained from the 
veil region; the reflectivity data demonstrate the reduced optical density of the dendritic network 
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as compared with adhesive regions and the topography data indicate that the veil is nearly 50 
nm above the adhesive regions. Finally, the data obtained from these dynamic filopodia suggest 
that retrograde flow is the dominate mechanism responsible for retraction although disassembly 
at the tip likely contributes [61-64]. 

Live cells can exhibit a range of motion, from extremely rapid to essentially stationary in 
the field-of-view. Therefore, the applicability of this method depends both upon the rate of 
change of the phase at a given pixel and the rate of change of the phase between adjacent 
pixels. As with other forms of interferometry, for large changes in optical path length, if the 
phase is under-sampled, it will not be possible to successfully unwrap the 2n phase changes 
should these occur. Because the laser feedback interferometer has been demonstrated to respond 
linearly and without distortion to dynamic changes in the optical path length until about 1 
MHz [30], it should be possible to use this approach to measure thermally and actively driven 
membrane fluctuations and dynamic changes in cell shape and organelles, an area where the 
application of RICM has been used extensively; it has been observed that these thermally driven 
membrane fluctuations at physiological temperatures are of order 10 nm rms with a wavelength 
<0.5 lim [65]. 

Alternatively, rather than image large changes in optical path length, psLFIM should be an 
ideal method for studying the earliest phases in the birth of integrin adhesions. In particular, 
we plan to apply the method in order to verify the predictions of a model that we have de- 
veloped that provides a mechanistic understanding of the processes that govern the earliest 
integrin adhesions [66]. In analogy with the cellular measurements presented in Ref. [67] and 
the measurements on vesicles decorated with receptor and/or repeller molecules [68-72], the 
sensitivity of psLFIM to static and dynamic changes in the actin network should permit its ap- 
plicability to the measurement of the nucleation, dissolution and merging of integrin adhesions 
at the ventral surface of cells on glass substrates. 

8. Conclusion 

We have coupled phase-shifted laser feedback interferometry with an inverted microscope and 
developed and verified the accuracy of an algorithm that addresses the inverse problem so that 
we can separate the contribution of the reflection from the coverslip-buffer interface from the 
field that has been elastically scattered by sub-cellular structures. Using this approach, we have 
reconstructed the topography and reflectivity of regions of the ventral membrane at focal ad- 
hesions. As expected, there is a strong correlation between the presence of paxillin and actin 
in stress fibers as identified with fluorescence imaging) and the interference images. Unlike the 
data obtained with RICM, which show dark regions at focal adhesions, using psLFIM we are 
able to demonstrate that these optically dense regions on the cytoplasmic face have high reflec- 
tivity and we are able to measure the topography of the regions independently of the reflectivity. 
The data acquired from fixed and live cells show the presence of a dense actin layer located «s 
100 nm above the coverslip interface. Additionally, the measured dynamics of filopodia and 
the lamella in a live cell supports retrograde flow as the dominate mechanism responsible for 
filopodia retraction. 
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